Macrospin approximation and quantum effects in models for magnetization reversal 
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The thermal activation of magnetization reversal in magnetic nanoparticles is controlled by the 
anisotropy-energy barrier. Using perturbation theory, exact diagonalization and stability analysis 
of the ferromagnetic spin-s Heisenberg model with coupling or single-site anisotropy, we study the 
effects of quantum fluctuations on the height of the energy barrier. Opposed to the classical case, 
there is no critical anisotropy strength discriminating between reversal via coherent rotation and via 
nucleation/domain-wall propagation. Quantum fluctuations are seen to lower the barrier depending 
on the anisotropy strength, dimensionality and system size and shape. In the weak-anisotropy limit, 
a macrospin model is shown to emerge as the effective low-energy theory where the microscopic spins 
are tightly aligned due to the ferromagnetic exchange. The calculation provides explicit expressions 
for the anisotropy parameter of the effective macrospin. We find a reduction of the anisotropy-energy 
barrier as compared to the classical high spin-s limit. 

PACS numbers: 75.60. Jk,75. 10. Jm 



I. INTRODUCTION 

Nanosystems of interacting magnetic moments have 
attracted much interest due to ongoing technological ad- 
vances in the field of magnetic data-storage devices. With 
decreasing system size, ferromagnetic particles on solid 
surfaces^ or ferromagnetic molecular magnets like Fes 
and Mni2j^ for example, are found in a single-domain 
state. The question of the stability of the ferromagnetic 
state against different kinds of thermal and quantum fluc- 
tuations is a crucial technological issue which at the same 
time provides serious challenges to a theoretical modeling 
and understanding. 

The magnetic properties of a ferromagnetic nanopar- 
ticle can be described by a Heisenberg model 



H, 



id 



to a model with only two degenerate ground states in 
which all spins arc pointing into the +z or — z direction, 
respectively. The anisotropic ferromagnetic Heisenberg 
model H = Hj + Hd in fact represents the simplest 
model to study magnetization reversal, i.e. the thermally 
activated transition between two ground states across 
an energy barrier induced by the anisotropy^ The rate 
of magnetization switching and thus the magnetic sta- 
bility obviously crucially depends on the height of the 
anisotropy-energy barrier AE. 

For weak anisotropy D <C J, it is self-evident to con- 
sider the spins as tightly coupled by the exchange inter- 
action and forming a huge macrospin, 



(3) 



the dynamics of which can be approximated by an effec- 
tive model 



with positive exchange coupling > between micro- 
scopic spins Si and Sj. The spins give rise to magnetic 
moments which are assumed to be localized at the sites 
i = l,...,L. The latter may constitute a ti-dimensional 
lattice of finite size with L sites in total and constant ex- 
change J = Jij between nearest neighbors i and j only. 
In a ground state of the model all spins are perfectly 
aligned, and the particle is ferromagnetic. However, this 
state is unstable against thermal fluctuations. Strictly 
speaking, for a system of finite size, the ferromagnetic 
state is destroyed at any finite temperature by fluctu- 
ations originating from a coupling of the system to an 
external heat bath. 

On the other hand, anisotropic contributions to the 
Hamiltonian, such as a single-site anisotropy of the form 



H,, 



Ht 



(2) 



which for D > distinguishes the z-axis as the easy axis, 
stabilize the ferromagnetic state. Inclusion of Hp leads 



^ macro ^ z 



const. 



with an anisotropy barrier 

A.E = -Dmacro'S' 



(4) 



(5) 



In this macrospin model, magnetization reversal takes 
place at zero temperature by suppression of the barrier 
due to an external magnetic field as described by the 
Stoner-Wohlfarth models At finite temperature T, re- 
versal may be caused by thermal activation. For low T, 
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FIG. 1: Two different models for magnetization reversal. A: 
Coherent rotation. B: Nucleation and domain-wall propaga- 
tion. Figure inspired by Ref. 
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the switching frequency is exponentially small and a fer- 
romagnetic state is stable over macroscopically relevant 
times while for high T, in the so-called superparamag- 
netic state, the magnetization vanishes on the time scale 
of the measurement. Classical theory^ for superparam- 
agnetic dynamics based on the Landau-Lifshitz equation^ 
and including a stochastic Langevin field to simulate a 
thermal bath leads to an Arrhenius law for the thermal 
activation rate T oc exp(— f3AE). This Neel-Brown law 
again emphasizes the role of the energy barrier AE for 
the magnetic stability of the particle. 

For stronger anisotropy, the macrospin approximation 
is no longer applicable as magnetization reversal prefer- 
ably takes place by nucleation and domain-wall propa- 
gation (see Fig. Q}. Here, classical many-spin models- 
are frequently used to study the reversal mechanism and 
in particular the transition from coherent rotation to 
nucleation and domain-wall propagation with increasing 
strength of the anisotropy. In particular, classical Monte- 
Carlo simulations are employed^ where Monte-Carlo up- 
date steps are related to physical time propagation. As 
a function of the model parameters, reversal mechanisms 
different from coherent rotation, can be identified and 
"phase boundaries" can be found in this way. The strong 
influence of the shape of the nanoparticles on the switch- 
ing rate, as observed experimentally,^^ can be explained 
by theoretical analysis of domain-wall propagation! 12 i 13 

The macrospin has a high spin quantum number S 
and may be treated to a good approximation classically. 
There are, however, different types of quantum effects to 
be considered: For example, quantum tunneling as a re- 
versal mechanism is competing with thermally induced 
reversal.— ~— In case of not too small systems, tunneling 
is significant only for systems extremely isolated from 
any dissipative environment at ultralow temperatures. 
For a macrospin coupled to a conduction band, Kondo 
screening is another issue. Since S > 1/2 there is an 
underscreened Kondo effect, the development of which is 
hampered by the anisotropy as has been shown by time- 
dependent numerical renormalization group.— As the 
quantum character of the macrospin is not accounted for 
by simple Langevin dynamics^ different approaches^— 
have been suggested to treat the macro- or many-spin dy- 
namics in contact with bosonic baths by means of quan- 
tum master-equation approaches and to determine, e.g., 
the blocking temperature. The importance of a many- 
spin model, for example, is demonstrated by studies of 
the probability density function of the macrospin ob- 
tained by replacing thermal with Markov processes^ 
Differences between the classical-spin limit and large 
quantum spins have been discussed^ 

Here we reconsider the question of competing reversal 
mechanisms by taking a quantum many-spin model as 
the starting point. The competition between coherent ro- 
tation and nucleation can be addressed in a less ambitious 
way by studying the static properties of an anisotropic 
quantum Hcisenbcrg model H = Hj + H ani . As the 
single-site anisotropy @, which is frequently studied in 



the classical context, cannot be used for the lowest spin 
s = 1/2 in a quantum model since trivially sf = const, 
we also consider a coupling anisotropy of the form 

#A = -g ^2&ijSizS jz , (6) 

ij 

i.e. .Hani = Hd or i? an i = Ha- This model will be inves- 
tigated by different complementary techniques including 
perturbation theory in -ff an ;, exact diagonalization and 
the Lanczos approach for systems with small size L and 
classical stability analysis in the high-spin limit. 

There are different goals of the present paper: As the 
status of the macrospin approximation is less clear in 
the quantum case, we aim at a strict derivation of the 
macrospin model in the weak-anisotropy limit. In this 
way it should be possible to relate the effective anisotropy 
strength of the macrospin to the parameters of the un- 
derlying many-spin model. The dependence on the spin 
quantum number s should exhibit quantum effects for 
small s and recover the classical result for s — > oo. Fur- 
thermore, beyond the weak-anisotropy limit, the break- 
down of the macrospin approximation and the competi- 
tion between the different reversal mechanisms, depend- 
ing on system size, coordination, dimensionality and on s, 
will be investigated. Our goal is to understand whether 
or not it makes a qualitative difference for the transi- 
tion from coherent rotation to nucleation with increasing 
anisotropy strength if starting from a classical or from 
a quantum spin model. Finally, the high-spin limit is 
addressed as a reference and for comparison with the 
quantum- mechanical calculations. 

The paper is organized as follows: In the next sec- 
tion [II] we briefly introduce the model, the basic concepts 
and notations by referring to exact-diagonalization re- 
sults. Sec. IIIII present the results of first-order pertur- 
bation theory for weak anisotropy and the derivation of 
the macrospin model. Its limitations are discussed in 
Sec. HVI on the basis of calculations employing the Lanc- 
zos technique. The dependence of the anisotropy-energy 
barrier on the spin-quantum number s and classical sta- 
bility analysis in the high-spin limit is addressed in Sec. 
IV1 Finally, Sec. I VII concludes the paper. 

II. EXACT-DIAGONALIZATION RESULTS FOR 
THE ANISOTROPY BARRIER 

Due to the SU(2) spin-rotation symmetry of the 
Heisenberg model, Eq. (fTJ , the total (or macro) spin, Eq. 
©, is a conserved quantity: [S, H] = 0. We consider the 
eigenstates of Hj: 

Hj\n,S,M) = E n (S)\n,S,M) . (7) 

The states are classified according to their total spin 
quantum number S — S m - m , ...,jf? max — l,S mgx and their 
total magnetic quantum number M = —S,...,S— 1,5. 
Here, Smax = Ls and S m - m = for s even or L even, 
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while iSmm = 1/2 otherwise. Throughout of the paper, 
we assume an even number of sites L for the sake of sim- 
plicity. Further, n labels the eigenstates in the invariant 
subspace with fixed S and M. For a given S, the SU(2) 
invariance implies the eigenenergies E = E n (S) to be 
25 + 1-fold degenerate and independent of M. 

For the present case of ferromagnetic exchange cou- 
pling Jij > 0, the total spin quantum number is maxi- 
mal, S = Smax = Ls, for a ground state. It is easy to 
see that the fully polarized state with all spins pointing 
into +z direction, \F) = \m\ = s) ® ■ ■ ■ ® \m L = s), is 
a ground state. Assuming all exchange interactions to 
be equal between nearest neighbors, Jij = J > 0, the 
ground-state energy is given by E (S max ) = ~(J/2)Ls 2 z 
where ~z = z i/ L is the average of the coordination 

numbers £j. 

The 2S , max + l = 2Ls+ 1-fold degeneracy of the ground- 
state energy will partly be lifted by the anisotropy if an i, 
Eq. (JU or Eq. flBJ). H ani breaks the SU(2) symmetry, 
[S 2 ,H an i] ^ 0, but preserves the rotational symmetry 
with respect to the easy axis (z axis): [S Zl H an {\ = 0. 
The anisotropic eigenstates |m, M) are therefore charac- 
terized by M and by an additional quantum number m 
labeling the states in the subspace with fixed M. The 
corresponding eigenenergies E m (M) will depend on M 
with 



E m (M)=E m (-M) , 



(8) 



resulting in a symmetric anisotropy barrier with a height 
AE given by the respective ground-state energies of the 
full Hamiltonian in the invariant subspaces with M = 



AE = E (M = 0) - E Q (M = S max ) > 



(9) 



For systems with a moderately large Hilbcrt space, i.e. 
for small s and small number of lattice sites L, the energy 
eigenvalue problem can be solved numerically. In case of 
the anisotropic problem Hj + -ff an i, conservation of S z 
can be exploited only. The dimension of the invariant 
subspace H(M) is given by^^ 



dimH(M) 



[(S max -M)/(2s+l)] 
k=0 

L-l+S max -M-{2s + l)k 
L- 1 



(10) 



where the binomial coefficients 



= for b > a, and 



where [x] denotes the greatest integer with [x] < x. For 
s = 1/2 full diagonalization is easily possible for systems 
with L = 14 spins. 

A simple example for a spin-s = 1/2 Hcisenberg chain 
with and without coupling anisotropy A on L = 10 
sites with open boundaries is given by Figs. ©, §S§ 
and (HJ). The perturbed (A > 0) ground-state energy 




FIG. 2: Energy eigenvalues of the isotropic s = 1/2 Heisen- 
berg model Eq. (TTJ on a linear chain of L = 10 sites and 
open boundaries as obtained numerically by exact diagonal- 
ization. Eigenvalues are classified according to the total mag- 
netic quantum number M = —Sm&x, Smax with S'max = 5. 
The nearest-neighbor exchange coupling J = 1 sets the energy 
scale. 



in the one-dimensional M = S mayi subspace is easily cal- 
culated as E (M = 5 max ) = -(J/2 + A/2)Ls 2 z. The 
A I = ±(S' max — 1) subspaces are reached by a single spin 
flip, i.e. excitation of a magnon. As there must be as 
many magnons as lattice sites, the dimension of the sub- 
spaces is given by L each. Their discrete energy spectrum 
can be seen in Figs. ((2J and (j4|) for A = and A = J, 
respectively. 

The low-energy sector of the isotropic spectrum for 
arbitrary M is shown in Fig. [3] Ground states have the 
maximal total spin S = 5 max . The ground-state energy 
is 2S' max + 1-fold degenerate. Excited states with M = 
±(5'ma X — 1), corresponding to magnon excitation, have 
total spin S = S max — 1. 

Fig. @] demonstrates the effect of an anisotropy term 
in the Hamiltonian. As the energy eigenstates are su- 



o -2.1- 









+ 


+ 

t 


+ 




+ 


+ 
+ 


+ 




+ + 

+s=t + 


+ 


+ 


+ 


+ 


+ 


+ + + +S=4 






+ 


+ 


+ 


+ +S=2 




+ 


+ 


+ 


+ 


+ + +S=3 


+ 


+ 


+ 


+ 


+ 


+ + + +S=4 


- + + 


+ 


+ 


+ 


+ 


+ + + + + S-5 

,1,1,1,1, 


6 -5 -4 


-3 


-2 


-1 





1 2 3 4 5 f 



FIG. 3: The same as Fig.[2]but on enlarged energy scale close 
to the ground-state energy. The total spin quantum numbers 
of the different spin multiplets are indicated. 
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FIG. 4: The same as Fig.[2]but with the coupling-anisotropy 
term, Eq. {BJ, included. A constant anisotropy Ay = A 
between nearest neighbors is assumed. Calculations have been 
performed for A — J. 



perpositions of states with different S (but the same M) 
in this case, their energy becomes M dependent. Conse- 
quently, an anisotropy barrier develops which, for a very 
strong anisotropy A = J, is large as compared to the 
finite-size gap between the ground states and the first 
excited states. Still, there is a residual degeneracy of 
the two ground states \F) and | — F) in the subspaccs 
A I = ±S max which implies that thermal fluctuations de- 
stroy ferromagnetic order. However, switching between 
\F) and | — F) now requires that thermal fluctuations 
must overcome the barrier. The anisotropy thus leads to 
a superparamagnetic stabilization of the magnetic state 
on a certain time scale as described by the Neel-Brown 
model.— i 



III. PERTURB ATIVE DERIVATION OF THE 
MACROSPIN MODEL 



In the weak- anisotropy limit, the different microscopic 
spins are tightly bound together by the ferromagnetic 
exchange coupling and form a huge macrospin which ro- 
tates from +z to — z direction, i.e. the microspins rotate 
coherently. To make this intuitive argument rigorous and 
quantitative, we derive the macrospin model as an effec- 
tive low-energy model by means of first-order perturba- 
tion theory in if an i in the following. This will result in 
a linear dependence of AE on the anisotropy strength. 
Non-trivial dependencies, however, may be expected with 
respect to L, s and the matrix J^- of exchange-coupling 
constants. 

First-order perturbation theory requires to compute 
the matrix element (S max , M'\si z Sj Z \S max , M) where the 
states |S, M) form an orthonormal common eigenbasis of 
S 2 and S z . Note that this includes the case of coupling 
(i =/= j) and single-site (i = j) single-site anisotropy. To 



find the action of Si Z Sj Z onto the state IS'max, M) we write 

V(Smax + M)\ c c _ M| 



ISmax.M) 



V / (^max-M)!(2S max )! 



Si™- M \F) (11) 



in terms of the fully polarized state with all spins point- 
ing in +z direction, \F) = |5 max ,M = S max ) = |mi = 
s) ■ ■ ■ \rriL = s). Here S_ = S x — iS y . Eq. ([TTj) is obtained 
by straightforward spin algebra. Using the commutator 

[Si z Sj z ., S_] = kSl_ (si—Sj z -\-Sj—Si z ) 



fc-1 



k 5% j S 

- D_ {Si-vj- 



-5^- 2 (.s l _ S ,_+ Sj _.s l _)(12) 



Si x ± isiy , we are left with matrix elements 



where s,± 

of the form (S max , M'\(D\F) where O contains microspin 
operators s... on the right-hand side acting on \F) as well 
as macrospin operators S... on the left-hand side which 
are considered to act on (S max ,M\. The effect of the 
former operators is trivial. The effect of the latter ones 
is obtained with the help of the relation 



7*| 



c,M) 



. V{s m ^ + M + ky.{s max -Aiy. 

V(S max + M)!(S max -M-fc)! 

X |S*max, M + k) • 



Using 



for M = S n 



(S max ,M| Si _|F) 



1 and 



2.s 



(S n 



Mis,: 



.\F) 



8s 2 



^4S max (2S 

max 

4s(2s- 1) 

max 



(13) 
(14) 

(i + 3) 

(i = 3) 
(15) 



for M = S max — 2, a straightforward calculation leads to 



(Smax, M\s iz S jz \S ma ,x,M') 



= S 



MM' 



s 2 + (2s-Sij) 



M - L z s 



2„2 1 



2L 2 s - L 



(16) 



This is a remarkably simple result which may be used 
to compute the anisotropy energy barrier Eq(M) in dif- 
ferent situations. To give an example, we consider the 
model H = Hj+H ari \ where a constant exchange interac- 
tion Jjj = J between nearest neighbors is assumed (and 
likewise for A y - in case of a coupling anisotropy). This 
implies that geometrical properties enter via the lattice 
topology only. For the case of the single-site anisotropy 
0,5^ 1/2, we get: 



E (M) 



1 , / 2s - 1 9 9 L(L - 1 
--JLs 2 z-D[ — M 2 + s 2 -±- '- 

2 \2Ls-l 2Ls-l 

0(D 2 ) , (17) 
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while for the coupling anisotropy (|6]) with non-zero and 
constant = A between nearest neighbors only, we 
find: 

E . ( M) . -ijiA-A^-^-y^) 

+ 0{A 2 ), (18) 

where z = z%/L is the average coordination number. 
This corresponds to a macrospin model Eq. (Q| with 

D_-£±D (19) 

for the single-site anisotropy, and 

~zs 

D macm = 777 T A ( 20 ) 

2.Ls — 1 

in case of the coupling anisotropy. The anisotropy en- 
ergy barrier is given by AE = D macm S 2 laK . Eqs. (JTSJ) 
and (|20[) provide an explicit and quantitative relation 
between the anisotropy parameter of the effective low- 
energy macrospin model and the microscopic model pa- 
rameters. 

Quite generally we note that this relation is non-trivial. 
In the limit L — > oo we find D ma cro oc l/L formally. 
However, for fixed fl or A, the macrospin approxima- 
tion becomes less accurate with increasing system size L: 
Eventually the finite-size gap, which controls the valid- 
ity of the perturbative treatment, is of the same order 
of magnitude as the anisotropy strength, i.e. the expan- 
sion parameter. The classical limit is reached for s — > oo. 
Here, L> macro = D/L or D macm = (z/2)A/L, respec- 
tively. The validity of the macrospin approximation for 
the classical model is discussed below. 



IV. LIMITATIONS OF THE MACROSPIN 
MODEL 

It is easily seen that the limits of weak and of strong 
anisotropy correspond to two qualitatively different ways 
of how to overcome the barrier. Consider the model sys- 
tem discussed above but with a large coupling anisotropy 
A> J. In the extreme case J = 0, we recover the Ising 
model, H = which exhibits a highly degenerate en- 
ergy spectrum. It is displayed in Fig. [5] Still there are 
two degenerate ground states in the M = ±S max sec- 
tors which are separated by an energy barrier. However, 
the barrier is entirely flat and thus qualitatively different 
from the quadratic trend given by the macrospin model. 
The tendency towards a flat barrier is already seen in 
Fig. |U for A = J. 

The energy eigenstates are constructed trivially in the 
Ising limit. In the ground state with M = S max the spins 
are aligned ferromagnetically, i.e. \S max , S max ) = | f>t 
,...,t>. Its energy is E (M = 5 max ) = -{A/2)Ls 2 z = 
—2.25 for the example shown in the figure. The ground 
state with M = S max — 1 is found by flipping a single 
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spin, i.e. \S niax , S max - 1) = | 4,t, — ,f). The correspond- 
ing excitation energy is AE = 2As 2 = 0.5, resulting in 
E (M = S max - 1) = -1.75. Shifting the "domain wall" 
to the right, i.e. \S max ,M) = | |, ...,|,t, -,t)> is possi- 
ble without further energy cost. Hence, nucleation and 
domain-wall propagation is the apparent mechanism for 
magnetization reversal in the Ising limit. 

The breakdown of the macrospin approximation hap- 
pens for still weak anisotropics but is gradual rather 
than abrupt. This is demonstrated by the exact- 
diagonalization results for the anisotropy-energy barrier 
Eq(M) shown in Fig. [5] For strong anisotropy A = J 
(corresponding to Fig. [4]), the macrospin approximation 
Eq. ([18]) completely fails and strongly overestimates the 
barrier height AE. With decreasing A, the macrospin 
approximation becomes more and more reliable. Still, 
at A — J/10 deviations from the exact-diagonalization 
results are visible on the scale of Fig. El If interpreted 
within the classical theory, this is still significant as the 
energy barrier affects the thermal activation rate expo- 
nentially strong. 

The fact that the breakdown of the macrospin model 
with increasing strength of the anisotropy is gradual 
rather than abrupt also implies a gradual change be- 
tween the corresponding reversal mechanisms. Within a 
quantum model, there are no well-defined phase bound- 
aries separating coherent rotation from nucleation and 
domain-wall propagation. This becomes obvious again 
in Fig. [7] where the anisotropy energy barrier AE is 
shown as a function of A. The macrospin and the exact- 
diagonalization results start to deviate from each other 
at arbitrarily weak A, i.e. already at second order in the 
anisotropy strength. 

Nevertheless, a typical anisotropy strength can be 
identified that marks the smooth crossover from coher- 
ent rotation to nucleation or the qualitative breakdown 
of the macrospin approximation. Fig. [7] shows that this 
crucially depends on the system size L. For L = 10 
spins (red lines) , the macrospin model is valid up to much 
stronger A as compared to the case with L = 24 spins 
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FIG. 6: Anisotropy barrier Eo(M) as obtained from the 
macrospin model (Eq. (|18p . lines) compared to the exact- 
diagonalization result (points) for three different anisotropy 
strengths as indicated. Results for an open chain of L = 10 
spins with s = 1/2 and nearest-neighbor exchange and cou- 
pling anisotropy. 



(blue lines). Generally, with increasing L, substantial 
deviations from the linear trend of AE(A) are found for 
weaker and weaker A. In the strong- anisotropy limit, 
AE(A) approaches a linear trend again which is given 
by the Ising limit (black line). Note that for the one- 
dimensional model considered here, the domain-wall en- 
ergy and thus A in the strong-anisotropy limit is indepen- 




FIG. 7: Height of the anisotropy energy barrier AE as a func- 
tion of the (coupling) anisotropy strength A for open chains 
of L = 10 and L = 24 spins with s = 1/2 as obtained by the 
macrospin model (blue lines, Eq. ([20])) and by the full model 
(red lines). Black line: AE as a function of A for the Ising 
limit with J = for both, L = 10 and L = 24. 
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FIG. 8: Height of the anisotropy energy barrier AE as a 
function of the (coupling) anisotropy strength A for systems 
of L = 24 spins with s = 1/2 as obtained by the macrospin 
model (Eq. ((20)) ) and by the full model (AE in the full model 
is always smaller). Results for different planar geometries 
as indicated. AE for the fully connected model, where the 
macrospin approximation is exact, is shown for comparison 
(black line). 



dent of the system size while in the macrospin or weak- 
anisotropy limit AE is roughly proportional to L (see Eq. 

as])). 

Calculations for L = 24 spins with s = 1/2 can no 
longer be carried out by full diagonalization. The present 
results have been obtained by employing the Lanczos 
algorithm ! 23 ' 24 This approximates the ground-state en- 
ergy Eq(M) of H in the invariant subspace H(M) by the 
ground-state energy of H in a Krylov space 

K n (M) = spanll^M^FIi.M),...,^"- 1 !^)} c H(M) 

(21) 

of dimension n <C dimH(M) that is spanned by 
the states \i, M),H\i, M), iT™ _1 |i, M) where \i,M) £ 
T~L(M) is an arbitrary initial state. Typically, n ~ 100 
Krylov-space dimensions arc fully sufficient to obtain an 
excellent approximation which, on the scale of the figures, 
cannot be distinguished from the exact result. 

The extent of applicability of the macrospin model 
strongly depends on the system's geometry. It turns out 
to become more and more reliable with increasing average 
coordination number z. This is demonstrated with Fig. [5] 
where again macrospin and cxact-diagonalization (Lanc- 
zos) results for the A dependence of the energy barrier are 
compared. Here, the system size is kept fixed to L = 24 
while the system geometry varies from a one-dimensional 
chain to a compact two-dimensional array, i.e. the aver- 
age coordination number z increases. The corresponding 
increasingly better agreement of the macrospin approxi- 
mation with the exact-diagonalization results is easily un- 



7 



derstood by consideration of the extreme case (see black 
line in Fig. [5]): For Zj = L — 1 = z, i.e. for the fully 
connected model with J^- = J and = A for arbitrary 
pairs i,j, the Hamiltonian reduces to 

H = - JS 2 - A S 2 Z + const (22) 

and thus 

Eo(M) = -A M 2 + const , (23) 

i.e. in this limit the macrospin approximation is trivially 
exact. 



V. CRITICAL ANISOTROPY STRENGTH IN 
THE CLASSICAL LIMIT 

The limit of the classical anisotropic Hcisenberg model 
should be recovered for spin quantum number s — > oo. 
Fig. [9] shows the anisotropy barrier AE as a function of 
the coupling anisotropy strength A as obtained from the 
Lanczos technique for a small chain of L — 4 sites. This 
permits a study of the spin model with quantum numbers 
up to s = 10. To compare the results for different s, we 
consider rescaled parameters, 

Jij H> Jij/s(s + 1) , 
Aij k> Ay/ s ( s + l) , (24) 

which corresponds to a rescaling Si i-> s,/ \/s(s + 1) of 
the spin variables. Disregarding quantum fluctuations, 
this amounts to a classical spin of unit length |s»| = 1 in 
the limit s — > oo. 

For s = 1/2 we again note the smooth crossover from 
the macrospin limit and coherent rotation at weak A to 
the Ising limit and nuclcation and domain-wall propaga- 
tion at strong A. In the quantum model the breakdown 
of the macrospin approximation starts immediately, i.e. 
the linear-in-A trend of the barrier is hardly visible as is 
apparent from the inset of Fig. [5J The crossover to the 
Ising limit takes place quickly. 

For s = 1, the range of anisotropy strengths where an 
almost linear trend is found, increases while the crossover 
to the Ising limit is seen to be delayed. With increasing 
spin quantum number s, these changes continue in a sys- 
tematic way and become more and more pronounced. 

From the result for s = 10, one can easily anticipate the 
s^oo limit: Here the anisotropy barrier is strictly linear 
in A for weak anisotropies up to a certain critical value 
A c /J = 1. At A c wc find a cusp, i.e. a discontinuous 
jump of d 2 AE/dA 2 . Beyond this point, a fairly broad 
crossover region follows until finally a linear-in-A trend 
is established again in the A — > oo limit. 

The strong-anisotropy limit is easily understood, as 
quantum fluctuations are suppressed anyway for all s. 
For weak anisotropy, on the other hand, the macrospin 
approximation is found to be exact in a finite range up 
to A = A c , i.e. perturbativc corrections of second and 



1 2 3 4 5 




A/J 



FIG. 9: A dependence of energy barrier AE for open 
chains with L = 4 spins and different spin quantum num- 
bers s as indicated. Exact results for the full model 
H.j + Ha after rescaling s; — > s,/ s(s + 1) for s — 
1/2,1,3/2,2,5/2,3,7/2,4,5,10. The inset shows the first 
derivative of AE with respect to A. 



higher order in A are decreasing with increasing s and 
finally vanish exactly in the classical limit s — > oo. For 
all A < A c , magnetization reversal takes place by a co- 
herent rotation of completely correlated spins in the clas- 
sical case. For any finite s, on the other hand, quantum 
fluctuations break up this perfect correlation. The mag- 
netic moments (sj) arc still aligned to the z axis during 
the reversal - we have (sj) oc e z since S z is a conserved 
quantity but the effect of fluctuations can be seen in 
the spin-spin correlation functions: In the M = sector, 
for example, we find (si z Sj Z ) = classically (s — > oo) and 
(sizSjz) = s 2 (l — 1/(1 — l/2Ls)) in the quantum case for 
sufficiently small A (see Eq. (fTl))) ). i.e. a small quantity of 
the order 0(L~ V ) which is independent of i and j. With 
increasing A at finite s, however, (si Z Sj Z ) gradually be- 
comes more and more site dependent. As a function of j 
and fixing i to one of the chain edges, (si Z Sj Z ) is monoton- 
ically decreasing in the M = sector. Spins at the chain 
edges are oriented antiferromagnetically: (si Z Sj Z ) < for 
2 = 1 and j = L. Hence, for any finite A in the quan- 
tum case there are already indications for domain-wall 
formation while in the classical limit there is a clear-cut 
distinction, given by a critical coupling A c , between the 
two reversal mechanisms. 

To understand the origin of the "phase transition" 
in the classical limit and to find a means to compute 
A c for arbitrary system size and geometry, we first 
performed numerical minimizations of the total energy 
E({si, sl}) of the classical model H = Hj + i? an i for 
both, -ffani = Ha and -ffani = Hp. Arbitrary spin con- 
figurations {si,...,s^} with = 1 and satisfying the 
constraint J2i=i s iz = M are considered. For M = this 
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FIG. 10: Sketch of the different magnetization-reversal 
mechanisms in the classical model for anisotropy strength D 
smaller (macrospin model, coherent rotation) or larger (nuclc- 
ation) than the critical value D c and for infinite D (nucleation, 
vanishing domain-wall width) as verified by calculations for 
L = 4 spins. 



provides the classical barrier height and the correspond- 
ing optimal spin configuration. Calculations have at first 
been done for different small system sizes L, typically 
L = 4, to reproduce the critical point anticipated from 
the results of Fig. [5] 

Wc found the optimal spin configuration with mini- 
mal total energy for a given M to be coplanar in all 
cases, i.e. there is a coordinate frame with Si y = for all 
spins. As dictated by the symmetry of the model, con- 
figurations that transform into each other by rotations 
around the z axis are degenerate. For a chain geometry, 
there is also degeneracy between configurations obtained 
by mirror transformations, Sj — > s^+i—i, and nucleation 
can thus start from both edges with equal probability. 
Furthermore, for M = wc have Sj. z = — S£_i + i j2 in the 
optimal configuration. The results can be summarized by 
the sketch of the optimal spin configurations at a given 
M in Fig. [TO] Results for coupling and for single-site 
anisotropy are found to be qualitatively the same. As 
is seen from Fig. I1Q[ the different mechanisms for mag- 
netization reversal in the two limits of weak and strong 
anisotropy are reproduced from the calculations for small 
L. 

More important, however, the physical origin of the 
critical point can be identified easily: At A = A c (or at 
D = D c , respectively), and approaching the critical point 
from the weak-anisotropy side, the optimal spin configu- 
ration exhibits an instability towards a non-collinear or- 
dering. As sketched in Fig. [TU1 the instability develops 
simultaneously for any M. 

This observation can be exploited to calculate the criti- 
cal anisotropy strength for arbitrary system size and sys- 
tem geometry in the following way: Using Si y = 0, we 
parameterize the spin variables as Si = (yT — mf, 0, m,) 
with — 1 < n%i < 1. This yields the energy functional in 
the form 

E({mi,...,m L }) = -g Aj-m^mj 



which comprises both, the coupling and the single-site 
anisotropy case. To find the critical point, it is sufficient 
to concentrate on M = and to compute the magnetic 
response of the system in the M = state with spin 
configuration mi = • • • = mj, = 0. Close to the criti- 
cal point, we can assume m., -C 1. Expanding the total 
energy in powers of m, then yields 



E{{mx, m L }) = const. + - y Xij™i™>i + ®{ m f 



with the inverse susceptibility matrix 



Xi 



(26) 
(27) 



and where we have defined J; = J^j Jij ■ We assume the 
model parameters as real and symmetric, = Jji and 

Aj = Av 

Note that approaching the isotropic limit with 
anisotropy strength D — > 0, the susceptibility matrix ex- 
hibits an eigenvalue approaching zero. This corresponds 
to an instability of the state towards ferromagnetic order- 
ing with all mi > and indicates the incipient degeneracy 
of states with different M. 

For finite D, this eigenvalue is negative and thus im- 
plies an indefinite susceptibility. This indicates that, by 
construction, x refers to the thermodynamically unstable 
excited state with M = and m t = 0. The local insta- 
bility of this state towards a non-collinear alignment of 
the magnetic moments at the critical anisotropy strength 
D c shows up as a zero of another eigenvalue of x- F° r 
a short chain of L = 4 spins with coupling anisotropy 
this happens exactly at A c = J, in agreement with the 
Lanczos results displayed in Fig. |H1 The eigenvector cor- 
responding to the zero eigenvalue describes the critical 
spin profile which agrees with the M = spin configura- 
tion for D > D c sketched in Fig. [TUJ 

Results for the critical anisotropy strength as a func- 
tion of the system size are shown in Fig. [TT] for both, 
coupling and single-site anisotropy and for systems of dif- 
ferent dimensions d. In all cases we considered systems 
with open boundaries. Generally, the anisotropy strength 
up to which the macrospin approximation is found to 
be valid, decreases rapidly with system size but depends 
only weakly on the type of the anisotropy term. There is 
also a weak dependence on the form of the particle as is 
obvious from the comparison of cubic arrays and spheres 
in d = 3 dimensions. 

For small nanosystems, the dimension and shape de- 
pendence of D c is somewhat irregular. For a larger num- 
ber of classical spins, on the other hand, the critical 
anisotropy strength for which the macrospin model can 
be applied, is the stronger the higher the dimension of the 
nanostructure. In fact, in the limit L — > oo, the different 
results are found to follow a simple power law, 



D C (L) oc L- 2/d = r 



(28) 



where d is the dimensionality, L the number of spins and 
/ = L 1 ^ the linear extension of the system. 
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FIG. 11: System-size dependence of the critical value for 
the anisotropy strength D c in units of J for one-dimensional 
open chains, for two-dimensional square arrays and for three- 
dimensional simple cubic arrays and spheres with open bound- 
aries. Blue lines: single-site anisotropy. Red lines: coupling 
anisotropy. 



This can easily be explained by comparing the en- 
ergy barrier for coherent rotation A_E corro t. with the bar- 
rier for nucleation and domain- wall propagation AE nuc i,. 
While the volume of the system oc l d is relevant for 
AE cor rot oc Dl d , the d — 1-dimensional domain wall of 
width a(D) oc \J J/ D 25 i 26 essentially determines AE oc 
Da(D)l d ~ 1 . The critical anisotropy is then obtained 
from D c l d oc D c y/ J/D c I , resulting in the scaling law 

(El. 



VI. CONCLUSIONS 

Anisotropy terms coupled to the ferromagnetic Heisen- 
berg model of many-spin nanosystems give rise to an 
energy barrier AE that crucially determines the ther- 
mal activation of magnetization reversal and thus the 
magnetic stability of nanosystems. Employing perturba- 
tion theory, exact diagonalization and Lanczos as well 
as classical total-energy minimization and stability anal- 
ysis, different quantum effects have been revealed that 
substantially affect the barrier height. 

Quantum effects are absent in the limit of strong 
anisotropy D. For D/J ^ oo the Ising model is ap- 
proached, i.e. a classical model where there are no quan- 
tum fluctuations. The magnetization-reversal mecha- 
nism in the Ising limit is given by nucleation at the 
edge or surface of the system followed by the propaga- 
tion of a domain wall with minimal (unit) width. With 
decreasing anisotropy strength, the domain-wall width 
increases. Finally, a critical value D = D c is reached 
where the reversal mechanism changes to a coherent ro- 
tation of the spins that are tightly aligned by the ex- 
change coupling J and form a huge macrospin. This 



critical anisotropy strength is well defined in the clas- 
sical limit s — > oo (with model parameters rescaled by 
the factor s(s + 1)) and can be computed for arbitrary 
dimension, system size and shape by an analysis of the 
local stability of the susceptibility matrix in the globally 
(thermodynamically) unstable excited state at the top 
of the barrier. For D < D c , i.e. for coherent rotation, 
the classical anisotropy-barrier height is simply given by 



AE C 
i.e.: 



E(mi 



mi = 0)—E(mi 



1 



AE, 



class 



LD 



(29) 



This is the same as the barrier of a classical macrospin 
of length S = Ls, composed of microscopic spins of 
unit length |s| = 1, as described by a macrospin model 
H = D macro S 2 with a barrier AE = D mliCIO L 2 if 
-Dmacro = D/L. Classically, and in the regime where 
coherent rotations takes place, the macrospin approxi- 
mation is exact. 

The first quantum effect manifests itself in the absence 
of the a well-defined critical anisotropy strength. As 
can be seen in the spin-spin correlation function (si Z Sj Z ), 
there is rather a smooth crossover from nucleation and 
domain-wall propagation to coherent rotation with de- 
creasing D. This crossover becomes less and less sharp 
with decreasing s, i.e. as the importance of quantum fluc- 
tuations increases. 

Second, for any finite anisotropy strength, the 
macrospin approximation is no longer correct in the 
quantum case. Only within perturbation theory to first 
order in D, the macrospin model emerges as the effective 
low-energy theory of the anisotropic quantum Heisenberg 
model. Comparing the true anisotropy-energy barrier 
AE of the many-spin system with the one of the cor- 
responding effective macrospin model, we find 



AE < AE„ 



(30) 



In the quantum case, there is virtually never a reversal 
by pure coherent rotation but always some "admixture" 
of nucleation and domain-wall propagation. The latter 
leads to the decrease of AE as compared to Ai? macro . 
This quantum effect becomes more and more pronounced 
with increasing anisotropy strength, with increasing sys- 
tem size, with decreasing dimensionality and with de- 
creasing average coordination number. 

Third, the quantum derivation of the macrospin model 
by means of first-order perturbation theory provides an 
explicit expression for the anisotropy parameter [see Eqs. 
(HHJ) and (|2"0|) ]. This leads to an energy barrier that is 
lower than the classical one: 



AE macYO <c AE c \ ass 



(31) 



Consider a single-site anisotropy case as an example: Af- 
ter the rescaling D — > D/s(s + 1) we have AS macro = 
[(2s - l)/(2£a - \)][L 2 s 2 /s{s + 1)]D. Therewith, the 
classical result (|2"9"|) is recovered in the high-spin limit: 
lim^oo A£ macro = Ai? c i ass . In the low-spin case, e.g. for 
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s = 1 , which is the extreme case for single-site anisotropy, 
however, we have AE macro = (L/4)D. This is smaller by 
a factor of 4 than the classical barrier AE c \ ass . Even in 
the limit of large systems L — > oo, there is substantial 
difference, 

1 2s — 1 

A^macro = 7T ~~rLD < LD = AE c \ ass , (32) 

2 s + 1 

compared to the classical limit which is approached 
rather slowly, i.e. AE macro / AE c \ ass = 1 - 3s/ s + 0(s~ 2 ). 
We conclude that, in addition to quantum many-spin ef- 



fects, there is a quantum correction of the macro-spin 
model itself which, for low spin quantum numbers, can 
be as important as the former one. 
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